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ABSTRACT 

Observationally confirming spatial homogeneity on sufficiently large cosmological scales is of 
importance to test one of the underpinning assumptions of cosmology, and is also imperative 
for correctly interpreting dark energy. A challenging aspect of this is that homogeneity must 
be probed inside our past lightcone, while observations take place on the lightcone. The star 
formation history (SFH) in the galaxy fossil record provides a novel way to do this. We calculate 
the SFH of stacked Luminous Red Galaxy (LRG) spectra obtained from the Sloan Digital Sky 
Survey. We divide the LRG sample into 12 equal area contiguous sky patches and 10 redshift 
slices (0.2 < z < 0.5), which correspond to 120 blocks of volume ^0.04 Gpc 3 . Using the SFH in 
a time period which samples the history of the Universe between look-back times 11.5 to 13.4 
Gyrs as a proxy for homogeneity, we calculate the posterior distribution for the excess large-scale 
variance due to inhomogeneity, and find that the most likely solution is no extra variance at all. 
At 95% credibility, there is no evidence of deviations larger than 5.8%. 

Subject headings: cosmology: theory, large-scale structure of universe, early universe 
1. Introduction 



The ACDM concordance model is extremely 
successful, as it can fit most cosmo logical observa- 
tions with just 6 free parameters ( Komatsu et al.l 
201 if ). Testing the assumptions that go into this 
model is vital, but it is often neglected. In par- 
ticular, the model rests on the assumption of spa- 
tial homogeneity and isotropy on sufficiently large 



scales (for a review see Clarkson fe Maartensl20ld 
Maartensl 1201 ll : IClarksonll2012l) . It is therefore 
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appropriate and timely to devise observational 
tests that allow us to probe the homogeneity and 
isotropy assumptions. We know the isotropy as- 
sumption is well supported by detailed observa- 
tions of the cosmic microwave background, which 
has shown that temperature variations are only 
one part in 10 5 across the sky. However, ho- 
mogeneity is much more difficult to probe. Ho- 
mogeneity is not established by observations of 
the CMB and the galaxy distribution: we can- 
not directly observe homogeneity, since we observe 
down our past light-cone, recording properties on 
2-spheres of constant redshift and not on spatial 
surfaces that intersect that lightcone. What these 
observations can directly probe is isotropy. In or- 
der to link isotropy to homogeneity, we have to 
assume the Copernican Principle, i.e. that we are 
not at a special position in the Universe. The 
Copernican Principle is not observationally based; 
it is an expression of the intrinsic limitation of ob- 
serving from one space-time location. 

The importance of testing the homogeneity as- 
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sumption has been highlighted by the develop- 
ment of inhomogeneous 'void' models which can 
potentially explain apparent acceleration without 
any exotic physics. By changing the mean den- 
sity and expansion rate radially away from us, 
observations such as SNIa can be accommodated 



have been developed to ex t ract this information 



without any dark ene r gy (s e e e.g.. iBiswas et al 
2010t iMarra fe Notarl fcoilt IClarksonl I2012L for 
reviews.). However, it is difficult to fit all ob- 
servations - in particular the combination of Hq 
and the CMB - without requiring significant inho- 
mogeneity or other departure s from the standard 
mode l at early times as well (INadathur fe Sa rkar 
20 lit iBull et al.1 l2012t IClarkson fc Regis! 1201 it 



de Putter et al 



20121) . 



This implies that tests for homogeneity must be 
made throughout the history of the universe. Con- 
sistency tests which could uncover deviations from 
homogeneity can be used to prob e consistency of 



obser vables on our past lightcone (jClarkson et al 
120081) . Testing for the transition to homogeneity in 
the galaxy distribution on the lightcone, while as- 
suming a Fri edmann background, is another con- 
sistency test ( Scrimgeour et al"]|2012l) . Probing in- 
side our past lightcone is harder, however, because 
we cannot observe it directly. One method is to 
use the Sunyaev-Zel'dovich effect to observe CMB 
anisotropies from distant clusters (lGoodmanlll995t 
Caldwell fe Stebbinsl l2008t lMaartensll201ll ). An- 



other is to probe the thermal his tory in widely- 
separ ated regions of the universe ( Bonnor fc Ellis! 
19861) . as it should of course be the same in the 



standard model. 

In this letter, we apply a new method of test- 
ing homogeneity in the interior of our past light- 
cone for the first time, by comparing the fossil 
record of galaxies at different redshifts at dif- 
ferent times along their past world-lines, thus 
accessing different patches of the Universe at 
the same cosmic time. A full proof of homo- 
geneity would entail establishing homogeneity of 
the metric tensor. Here we apply a consistency 
test to check for violations of homogeneity, us- 
ing the star for mation rate as a pr obe, follow- 
ing the idea in iHeavens et al.l ( 20 111 ) . The fossil 
record, or the star formation history (hereafter 
SFH), can be obtained by analysing the shape of 
the galaxy spectrum, which encodes information 
about the histories of the component stellar pop- 
ulations, dust, and star formation. Various tools 



(e.g.. IHeavens et a 
12003 : lOcvirk et alT 



2000; 



2006; Vincoletto et al 



Cid Fernandes et al.. 

mm. 



of which we use the V Ersatile SPectral Analysis 1 
(hereafter VESPA, see iToieiro et all 120071 1200a 
for more details). These approaches rely on the 
assumption that the evolution of the stellar pop- 
ulations is well understood and that the current 
modelling of stellar population is accurate. We 
use VESPA to obtain the SFH within the time 
bin 11.5 to 13.4 Gyrs of stacked Luminous Red 
Galaxy (LRG) spectra located at different posi- 
tions on the sky and at different redshifts. We 
compare the histories of different patches of the 
Universe, using the local star formation rate as a 
proxy for homogeneity. 

This letter is organised as follows; in $2] we 
briefly describe the applicability of the VESPA 
routine as a test of homogeneity, and then describe 
the data, sstar formationimulated data and our 
method in Sj3] We present the results in S|4] and 
conclude in <|5] 

To calculate distances, and to map from red- 
shift to time, we ass ume a fiducial flat AC DM with 
best fit WMAP7 feomatsu et al.l l201lh cosmo- 



logical parameter values ^m, ^6> n s , Hq = 
0.729, 0.271, 0.045, 0.809, 0.966, 70.3 kms^Mpc -1 ). 
Since we are looking for deviations from homo- 
geneity, it is conservative to assume this rela- 
tion, whic h may be different in inhomogeneous 
universes ([Heavens et al.l 1201 II ). Any viable dark 



energy or modified gravity model will have a back- 
ground redshift-time relation that is close to the 
concordance model's. 

2. VESPA and homogeneity 

An illustrative diagram of our method is shown 
in FigHJ Here, for illustration only, we assume 
that LRGs form at a similar cosmic time, and have 
similar SFHs, which we illustrate by the galax- 
ies changing colour. VESPA recovers the SFH for 
each galaxy along its own world line (vertical lines 
in FigJT]), allowing us to compare the SFH at dif- 
ferent distances but at the same cosmic time, e.g. 
at positions 1,2,.., 5. In practice there is scatter in 
the SFH, due to sample variance on small scales 
and measurement error. We will consider these 



1 http: / /www-wfau. roe. ac.uk/vespa/ 
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later, and seek additional variance from large-scale 
inhomogeneity. 



3. Data and Method 

Data: All of the galaxies used in this study 
were drawn from Sloan Digital Sky Survey (see 



were drawn rrom Sloan Digital 5kv aurvey (see 
York et all l2000t iGunn et all 1200a ISmith et all 
20021 and references therein) Data Release 7 
(jAbazaiian et al.ll2009l hereafter SDSS DR7). We 
use 8.5 x 10 4 galaxies between the redshift range 
of 0.25 < z < 0.55, selected to b e Luminous Red 



Galaxies ([Eisenstein et al.l 1200 ll hereafter LRG) 



drawn from the VESPA database. 

We divide the SDSS survey footp rint into 12 
equa l area sky patches using HealPbQ (jGorski et al 



2005), and N z — 10 redshift slices, whose widths 
are shown in Table [TJ together with the total 
number of galaxies, and the approximate volume 
of the SDSS survey in each redshift slice. We 
hereafter refer to the galaxies in each sky patch 
at each redshift slice as a 'block' (B) of galaxies. 
We randomly select galaxies in each block into 
sub-samples of approximately 200 galaxies and 
stack the SDSS galaxy spectra for all galaxies in 
ea ch sub-sample follo wing the method presented 



in lToieiro et al.l (|201ll ). Stacking the LRG spectra 
(as opposed to averaging the SFHs of individual 
galaxies in a block) allows us to recover the aver- 
age star formation history of a block with higher 
resolution in lookba ck time (see discussi on in Sec- 



tions 3.2 and 3.3 of iToieiro et al.ll201ll ). We use 



VESPA to interpret the stacked spectrum in terms 
of a star formation and enforce VESPA to recover 
measurements in 16 time bins, r'. The time bins 
are in the rest-frame of the stacked spectra, or 
alternatively the rest-frame of the galaxy block 
Tb, and we refer to these quantities as the Star 
Formation Histories SFH(Tb,t'). Additionally, 
we enforce VESPA to only allow star formation 
in bins whose starting times are after the start 
of the Universe, calculated assuming our fiducial 
cosmology. 

Methodology: We next add the age of the 

Universe, calculated using our fiducial cosmol- 
ogy, at the average redshift of the galaxy block 
Tb, to the ages of the recovered VESPA bins for 
the stacked spectra. We map the values of the 




Big bang 



Fig. 1. — An illustration of the test of homogene- 
ity. We assume that LRGs form at a similar cos- 
mic time, and have similar stellar formation rate 
histories (SFH), which we illustrate by the galax- 
ies changing colour. VESPA recovers the SFH for 
each galaxy in the galaxies rest-frame. We re-bin 
the SFH to the common lookback time, and com- 
pare the local star formation rate, for example, 
locations 1,2, ..,5 to probe homogeneity. Galaxy 
world lines are shown in comoving coordinates. 



Redshift ID 


Range 


Ngals 


Total Vol Gpc a 


1 


0.200 < z < 0.279 


7874 


0.90 


2 


0.280 < z < 0.308 


9352 


0.46 


3 


0.309 < z < 0.327 


8532 


0.34 


4 


0.328 < z < 0.342 


8594 


0.29 


5 


0.343 < z < 0.359 


9181 


0.36 


6 


0.360 < z < 0.376 


8202 


0.39 


7 


0.377 < z < 0.398 


8754 


0.55 


8 


0.399 < z < 0.424 


8277 


0.71 


9 


0.425 < z < 0.457 


8272 


1.00 


10 


0.458 < z < 0.537 


8065 


2.91 



Table 1: The redshift identifier and range of the redshift 
slices, the number of SDSS LRGs within each slice, and the 
approximate total volume in Gpc 3 contained by the redshift 
slice. 



2 http: / /hcalpix.jpl. nasa.gov 
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SFH(Tb,t') to a common frame SFH(0,t) with 
bins denoted by r, the lookback time with respect 
to the current epoch, and have chosen the lowest 
bin to be at t = 0. When we map the VESPA 
time bins r' to the common time bin r we choose 
to maintain the bin widths, to avoid over-binning 
the data. Fig. [2] shows the start and end times 
of the rest-frame r' VESPA bins in lookback time 
(horizontal axis), for each redshift slice (vertical 
axis). The continuous solid black lines show the 
locations of the common frame r bins with start 
times greater than 2 Gyrs. LRGs form most of 
their stars very early on, and so we concentrate 
our test of homogeneity within the greatest time 
bin t = 15, which corresponds to a look-back time 
between 11.5 and 13.4 Gyrs. The reason for this 
is that all stacks have considerably (almost two 
orders of magnitude) more star formation in this 
time bin than others, and have smaller fractional 
errors (see below for the error assignment). The 
distribution of a large number of random variables 
(recall each stack has ~ 200 LRGs spectra) can 
be modelled as a Gaussian, following the central 
limit theorem. For lower r < 15 bins some stacks 
have zero estimated star formation and larger frac- 
tional errors, skewing the distribution, which can 
no longer be modelled by a Gaussian, so we are 
unable to define a robust likelihood function. 

The redshift slices of the blocks (recall that a 
block is a redshift slice/sky patch) are chosen to 
contain N s > 3 stacked spectra, which are con- 
structed from the (N s ) galaxy sub-samples within 
the block. For each block B, we calculate the av- 
erage Ab, and estimate the standard deviation of 
the block SFH, as, from the sub-samples. We 
determine the mean value \x of Ab, and further 
calculate the average value of Ab for all blocks at 
fixed redshift z, A z and the standard deviation of 
A z across the N z = 10 redshift slices, which we 
denote as a z = a[A z ). 

The dispersion a z is scatter arising from 
the re-binning of solutions SFHb^Tb, t 1 ) of 
blocks at different redshifts to the common frame 
SFHs i(0,r). Note that block-to-block inhomo- 
geneity would contribute to this, but only at the 
level of 1/12 of the variance, so it will affect our 
conclusions on the rms inhomogeneity by only 4%. 
We will, however, make this correction. 

We compute the Student i-distribution ts, for 



all blocks: 



f, = 



At 



(1) 



which determines the number of 'combined error' 
a = y/ 1 o 2 B + cr z , or the departure the measurement 
Ab, is from the mean or notional value for the 
entire sample at each time r. 

For illustration we compare the probability den- 
sity function obtained by the above analysis, with 
the theoretical probability density function f(t) of 
the t-distribution with r/ = 10 x 12 degrees of free- 
dom, which has the analytic form given by 



/(*) = 



1 



^5(1/2,77/2) 



1 + - 



-(*J+l)/2 



(2) 



where B(l/2,7j/2) is the Beta function. We see 
in Fig. 3 that the distribution of t s follows the 
expected distribution reasonably well. The grey 
shaded area shows the 95% range for t s statistics 
from 4000 Gaussian Random Samples of the SFH 
and errors. 

We now more formally model the data as hav- 
ing a gaussian distribution but with the possibil- 
ity of an extra fractional variance V arising from 
inhomogeneity. i.e. we assume homogeneity and 
check for consistency using the likelihood of the 
data given by 



P B (Y): 



1 



2na\ 



■ exp 



(A B - m) /24 . (3) 



K M 2 . (4) 



If we assume a uniform prior for V, then P(V) = 
HbPb(V) is the posterior for V given the entire 
block dataset. As a check, we show P(V) in Fig. 4 
for simulated datasets (sub-blocks) with variance 
N s a B , for different star formation histories: a con- 
tinuous SFH, a gaussian SFH with mean 10 and 
standard deviation V2 Gyr, an exponential SFH 
with a scale length of 0.5 Gyr, and a SFH equal 
to the mean of the data. SFHs are re-binned to 
the common frame. We see in all cases that the 
posterior is correctly maximised at zero, and an 
upper limit dependent on the SFH. 

4. Results 

In Fig 0] we show with the solid line the poste- 
rior distribution for the additional fractional vari- 
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Fig. 2. — The start and end positions of the rest- 
frame (r' in the text) VESPA bins in lookback 
time, for each redshift slice. The continuous solid 
lines black lines show the locations of the common 
frame (r) bins. We concentrate our test of homo- 
geneity within the hnal bin between 11.5 and 13.4 
Gyrs. 




t-distribution 

Fig. 3.— The Student-i distribution for the SDSS 
SFH data, assuming that the variance is due to 
a combination of small-scale (sub-block) sample 
variance, measurement error and scatter due to re- 
binning from the galaxy rest-frame to the present 
epoch. The vertical axis shows normalized fre- 
quency. The grey region give the 95% spread of 
Gaussian Random Samples from the data and er- 
rors. We over-plot the theoretical ^-distribution 
(dashed line). 



T 
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Fig. 4. — The ensemble probability that the dis- 
persion seen in the values of Ag are drawn from 
a Gaussian distribution around /Lt, as a function of 
an additional error component V, which is scaled 
by /j, 2 . The colored dashed lines show the proba- 
bility for each simulated SFH and the black solid 
line shows the data. 

ance V in the SFH of the blocks. The most proba- 
ble variance due to inhomogeneity is zero, and 95% 
of the posterior probability lies within V < 0.0032. 
Hence the 95% credibility interval for the addi- 
tional fractional inhomogeneity r.m.s., \fV (as- 
suming a uniform prior on V), is 5.6%, or 5.8% 
if we include a correction for the rebinning. The 
colored lines show the different sets of simulated 
data. We see that only the Uniform SFH has a 
peak which is not at V = 0, however it is con- 
sistent with at the < 95% confidence level. We 
note that by artificially reducing the additional 
variance on each block ere, the peak moves closer 
to V = 0. 

5. Conclusions 

Modern cosmology is built upon the assump- 
tion of homogeneity which is inferred through the 
observation of isotropy (e.g. by the Cosmic Mi- 
crowave Background radiation) and the Coperni- 
can principle, stating that we do not occupy a pre- 
ferred location. 

Deviations from homogeneity, in particular, 
an inhomogeneous background, for example LTB 
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models in which massive void exist, can poten- 
tially explain the dimming of distant supernovae 
without invoking dark energy. Testing homogene- 
ity is therefore an active area of research, and 
many tests have been devised, e.g., kinematic SZ 
effect. 

In this letter we have per formed a new obser - 
vational test of homogeneity ([Heavens et al]|201lf ) 
by examining the estimated Star Formation Histo- 
ries in old stars from s tacked spectra of SPSS Lu- 



200lh 



us- 



minous Red Galaxies (jEisenstein et al 
ing VESPA. The data are blocks in 10 redshift 
intervals 0.025 < z < 0.55, with 12 equal-area an- 
gular bins. 

We estimate the sample variance and measure- 
ment error arising from small-scale (sub-block) 
variations by computing the error on the mean of 
the sub-blocks. Additionally, we include the scat- 
ter arising from re-binning to the present-day look- 
back time, and then perform a Bayesian analysis of 
any additional variance which may exist on large 
scales. Our test assumes homogeneity and checks 
for consistency and we find no evidence for extra 
variance, and a 95% upper limit to the credibility 
interval of a fractional variation of 5.8% in SFH 
between 11.5 and 13.4 Gyrs. The typical block 
size is about 0.04 Gpc 3 . 

The main uncertainty is in the stellar pop- 
ulations models employed by VESPA. However, 
this result can be easily extended and improved 
upon with future spe ctroscopic surveys e.g., BOSS 



( Dawson et al. 20121 ). and as our knowledge of 
Stellar Population Models increases. Although 
this is not a complete test of homogeneity, which 
would require investigation of the metric tensor it- 
self, this limit on homogeneity is the first to come 
from within the past light cone, rather than being 
restricted to our past light cone. As such it is gen- 
uinely testing homogeneity rather than isotropy. 
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